Branching fraction of the isospin violating process <p — > cj7t' 
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We examine the parametrization of the e + e~ — > lutt cross section in the vicinity of the <f> reso- 
nance and the extraction of the branching fraction of the isospin violating process 4> — > < ^ 7r ° from 
experimental data. We found that there are two possible solutions of the branching fraction, one is 
4 x 10 -5 , and the other is 7 x 10 -3 . The latter is two orders of magnitude higher than the former, 
which is the commonly accepted one. 



PACS numbers: 14.40.Be, 13.25.Gv, 13. 66. Be 
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I. INTRODUCTION 

As has been pointed out in a recent study [l[ , there are 
many cases where multiple solutions are found in fitting 
one dimensional distribution with the coherent sum of 
several amplitudes and free relative phase between them. 
The fit to the e + e~ — > urn cross sections in the vicinity of 
the </> resonance was shown as an example of the existence 
of the two solutions and how large the difference could 
be between them. 

However, Ref. [l[ found these two solutions only 
through a fit to the experimental data, which causes sus- 
picion that the two solutions may due to the statistical 
fluctuation or other reasons associated with the data han- 
dling or fitting procedure or something else. In this brief 
report, we show mathematically that two solutions exist 
in the parametrization of the cross section used in the 
original publication Q; and the second solution can be 
obtained analytically from the solution reported in the 
literature without doing fit to the experimental data. 



II. DERIVE OF THE SECOND SOLUTION 



If we write 



The cross section of e + e — > 



center-of-mass energy, yfs : is parameterized as 



r(Vs) = CTnriVs) 



1 - z 



(1) 



in Ref. Q, where cr nr (y/s) = a + cr'{y/s — M^,) is 
the bare cross section for the non-resonant process, pa- 
rameterized as a linear function of y/s; M,p, T^, and 
— Mj — s — iAf^Tcf, are the mass, the width, and the 
inverse propagator of the <fi meson, respectively. Here Z 
is a complex number which depicts the interference ef- 
fect. Conventionally, the real and imaginary parts of Z 
are denoted as $t(Z) and ^s(Z), respectively. 
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G(s,Z) = 1-Z 



(2) 



then in the complex-parameter space (denoted by a com- 
plex number Z'), we want to figure out all possible pa- 
rameters which can satisfy the following relation 



\G(s,Z)\ 2 = \G(s,Z')r 



(3) 



Notice that the above relation is to be true for any s, it 
should be true for some special values of s. If we firstly 
take a special value of s which satisfies Mj — s — 0, then 
we obtain 



II -iZ\ 2 = \l-iZ'\ 2 



\Z\ Z + 23?(Z) = \Z'\ Z + 2<3(Z') 



(4) 



(5) 



Secondly, we take another special value of s which satis- 
fies M? — s = MfkTfk, we obtain 



1 + i 
1 —Z 



1 + i , 
1 —Z' 



(6) 



as a function of the or 



\Z\ 2 - 25R(Z) + 23(Z) = \Z'\ 2 - 2U{Z') + 2%Z') . (7) 

Subtraction of Eq. §5§ from Eq. ([7]) yields 

R(Z') = . (8) 

With this equality, Eq. ([5]) is recast as 

[1 + 3(Z')] 2 = [1 + 3(Z)] 2 , (9) 

by virtue of which one gets either ^s(Z') = ^s(Z) or 
Q(Z') = —2 — ^s(Z). As a summary, we have two sets of 
solutions: 



and 



K(Z') = R(Z) 



^t(Z') = , 
9f(Z') = -2-9f(Z) 



(10) 



(11) 
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It is readily to check that the above two sets of solutions 
are true for the relation for any s. Obviously, the first 
set of solution is trivial which can be expected intuitively. 
However, the second set of solution is fairly interesting 
which is firstly obtained analytically. More interesting 
thing is, according to the Eqs. (fTO]) and (fTT|) . the second 
set of solution can be obtained from the first one. Both 
solutions describe the experimental data identically well 
and one can not distinguish them purely from the ex- 
perimental data. Therefore we conclude that if the cross 
section of e + e~ — > ljtt as a function of the center-of- 
mass energy is parameterized as Eq. ([!}, there must be 
two sets of solutions of the interference parameter Z. 

One remark on our mathematical analysis. More gen- 
erally we write G(s, Z) in the form 



G(s,Z) = 1 + ZF{s) 



(12) 



with F(s) being a complex function depending on s. If 
two special values of s are taken, from relation ^ , we ob- 
tain two quadratic equations with two unknowns. Gener- 
ally speaking, there should be four sets of solutions. The 
existence of the exact two sets of solutions of our exam- 
ple indicates that the number of sets of solutions depends 
strongly on the form of F(s). 



III. EXPERIMENTAL CONFIRMATION 

From Ref. @, 3?(Z) = 0.106 and = -0.103 are 

acquired from a fit to the experimental data in the uj — > 
7r + 7T~7r° decay mode. In the light of Eq. (fTTj). the second 
set of solution can acquired immediately, i.e. %t(Z) = 
0.106 and = -1.897. 

It is interesting to compare these results with those 
obtained from a fit to the experimental data [ij. From 
Table 2 of Ref. |l| , we found that the real parts of the two 
solutions are identical, while the imaginary parts sum up 
to a number slightly different from —2. A check of the 
fit results showed in Ref. [l[ indicates that the —1.90 is 



indeed from a rounding of —1.897. Keeping one more 
digit, we find the sum of the imaginary parts is exactly 
—2. Both the real parts and the imaginary parts agree 
perfectly between the fit results and the analytical eval- 
uation. 



IV. SUMMARY AND DISCUSSIONS 



We showed above that there must be two solutions in 
extracting the branching fraction of — > 



parametrization of the e + e 



• cl>7t with the 
cross section around 
the 4> resonance in Ref. Q . While the first solution corre- 
sponds to B((f> — > W7T°) = 4 x 10~ 5 as reported in Ref. 0], 
the second solution would be B(<p — > ojit ) = 7 x 10~ 3 , 
which is two orders of magnitude higher than the first 
one. 

One may need to check whether the parametrization 
of the cross section is meaningful or if there are further 
constraints to the parametrization or the parameters, in 
order to pick out the physics solution from the two-fold 
ambiguities. 

It is worth pointing out that 4> —> is an isospin vi- 
olating process and thus should be small. The branching 
fraction reported in Ref. 2| is already large compared to 
theoretical calculations [3|. However, if the physics is the 
second solution showed above, we would find theoretical 
calculations arc too low. 
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